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by Richard J. Priem 
Lewis Research Center 

SUMMARY 

The relative importance of various steps in the combustion process on combustion 
instability was compared theoretically by assuming that individual steps controlled the 
growth of a disturbance. The combustion steps involved the time dependence of (1) pro- 
pellant injection, (2) accumulation of unburned propellant, (3) propellant atomization, 

(4) propellant vaporization, and (5) chemical reaction. The calculation gave stability 
limits (minimum disturbance that will grow to instability) with various steps considered 
to control the combustion rate. 

The calculations show that chemical reaction is the most sensitive mechanism, but 
if accumulation of unburned propellants is included, the chemical reaction step is likely 
to be important only for low conversion rates that do not occur in most rocket combus- 
tors. With realistic combustion conditions the physical processes of vaporization and 
atomization appear to be the most important in determining instability limits, the limits 
for each step depending on the accumulation of unburned propellant and vortex velocity. 
Experimental instability studies with an 8-foot torus combustor, which is similar in 
geometry to that used for the theoretical calculations, confirmed the calculated results. 
Experiments also confirmed qualitatively the fluctuations that occur during instability. 

INTRODUCTION 

Combustion of propellants in a liquid rocket combustion chamber involves the re- 
action of a reducing agent (fuel) and an oxidizing agent to produce reaction products with 
the liberation of heat. This seemingly simple statement of a reaction between the reduc- 
ing and oxidizing agents can become very complicated when the mechanism or steps are 
considered that occur between the injection of liquid propellants and the time or position 
that the final products are formed. Occurring simultaneously with the combustion are 
fluctuations in the chamber pressure, which, under certain conditions, are seen to 
amplify into periodic oscillations driven by the combustion process. These periodic 
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pressure oscillations result in engine failure due to ex- 
cessive heat transfer to the walls, which eventually 
leads to thermal failure, or vibrational interference with 
controls and other accessories that leads to mechanical 
failure. Thus, it is important to understand the role of 
the combustion process in the growth of these oscilla- 
tions. 


Figure 1. - Combustion section considered for Previous studies have shown how various parts of 

various instability models. , . , , . , , , , , 

the complicated combustion process could support the 
growth of a pressure disturbance. (See bibiography in ref. 1. ) These studies have 
usually resulted in qualitative descriptions of the combustion process with pressure 
oscillations or a theory for a particular step that could sustain an oscillation. In each 
study the assumptions, combustor design, model, etc. , are different; thus, it is im- 
possible to evaluate the importance of a particular step in the combustion process in sup- 
porting the growth of a pressure disturbance. Therefore, a theoretical study was con- 
ducted, similar to that reported in reference 2 which used a toroidal ring to represent 
the combustor, in which the various combustion steps were individually assumed to con- 
trol the growth of a disturbance. It was determined which step could amplify the small- 
est disturbance, so that this process could be further studied under more realistic and 
complicated conditions. To determine the influence of a vortex in the combustor, which 
may be produced by viscous forces with a transverse oscillation (ref. 3), calculations 
were also made with a constant vortex velocity superimposed on the oscillations. As a 
check on the theory, the calculated stability limits were compared to the experimental 
limits obtained with several injectors in a torus combustor. 


THEORY 

General Equations 

The combustor geometry and transport equations used in the analytical portion of 
this investigation are identical to those in reference 2. The combustor was an annular 
section with a very small thickness AR and length Az, as shown in figure 1. This 
geometry limits the oscillation to a transverse wave. The propellants are assumed to 
be uniformly injected and to burn at some rate which depends on their position within the 
combustor. The burning is accompanied by gas and liquid flow in the axial direction. 

As the propellants burn, random disturbances occur within the combustor. If the ampli- 
tude of a disturbance is large enough, it may develop into a wave in the combustor. The 
object of this study was to determine analytically, for various combustion steps, the 
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minimum amplitude of a pressure disturbance in the annulus that is needed to develop 
a wave. This calculation was accomplished by solving the transport equations that 
follow, which are derived in reference 2. Cylindrical coordinates are used in which the 
z -direction corresponds to the axial direction in the combustor. 
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(Symbols are defined in appendix A. ) The derivatives in the axial (z) direction were 
determined with the assumption that the total mass, momentum, and energy in the 
annulus were constant. Furthermore, it was assumed that these derivatives were in- 
dependent of e. These assumptions result in the following equations (derived in ref. 2), 
which permit the evaluation of the derivatives with respect to z: 
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Combustion Models for Burning Rate 

In addition to the transport equations, an expression for the instantaneous local 
burning rate is required. A schematic model for the entire combustion process is 
illustrated in figure 2. Propellants are assumed to enter the chamber at a varying rate 
oh dependent on the injector pressure drop P T - P c , which assumes no inertial or 
capacitance factors in the feed system. (The derivations of expressions for various 
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rates and masses are given in appendix B. ) Liquid is atomized at a rate dependent 

cl 

on th$ amount of unatomized liquid M . Similarly, the propellants vaporize at a rate 

cl 

o> v depending on the amount atomized but not vaporized M and then react chemically 
at a rate u> r depending on the amount of each propellant vaporized but unreacted M r - 
In figure 2 the arrows represent propellants changing from one state to another, and the 
crosshatched areas indicate propellants in a given state. The steps between propellant 
vaporization and reaction must include mixing, but since there are no known expressions 
to determine the local instantaneous mixing rate, it was not included in the model. 

To determine the relative importance of the various combustion steps in driving a 
pressure oscillation, several models were used. Each model assumed that one of the 
mechanisms controlled the local instantaneous combustion rate and that the other mecha- 
nisms exerted negligible effect on the overall rate; subsequent mechanisms were con- 
sidered very fast so that they would follow at the same rate. Preparatory mechanisms 
were considered fast enough to follow the injection rate and yet be insensitive to local 
fluctuations. With each model an adiabatic pressure disturbance was applied to the sys- 
tem, and its growth, or decay, was observed. This disturbance was given by 
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Chemical Reaction Model 


For this model it was assumed that only the chemical reaction rate varied during in- 
stability. The burning rate expression was then 
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Vaporization Model 

For this model it was assumed that the vaporization rates for oxidizer and fuel were 
equal and were the slowest compared to those of the other rates. It was also assumed 
that the atomization and chemical reaction mechanisms were very rapid so that they 
followed the injection and vaporization rates, respectively, or 
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The burning rate expression then becomes 
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Atomization Model 

It was assumed that the atomization rates for oxidizer and fuel were equal and were 
the slowest for this model and that the vaporization and chemical reaction followed them, 
or 
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The burning rate was then obtained from 


w' = or 


a, ox 


co T , = M» (p’) 5//12 
a, f a v ’ 


otO. 625 


1 + 


Av'/ 


M a,ox = M a,f=^^ V 

a, o 


f («• - <^)dt' 

Jo 


(15) 


(16) 


li 




. .x 

Burning rate parameter, $ 

Figure 3. - Instability boundaries for chemical 
reaction model with finite concentrations. 



THEORETICAL RESULTS 

The results of the numerical calculations were pressure, particle velocity, density, 
temperature, burning rate, and concentration of mass at each of 20 positions in the 
annulus as functions of time. Typical plots of pressure and velocity can be found in 
reference 2. From these types of curves it was possible to determine the amplitude of 
the minimum disturbance that would increase with time. This value of the minimum dis- 
turbance was used to determine the boundaries between stable and unstable combustion. 
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All calculations presented herein were for the following conditions of v^, AP', 
Av’, and y in equations (1) to (17): 


Viscous -dissipation parameter,^ 0.3x10”^ 

Nondimensional axial gas velocity, v* 0. 05 

Nondimensional steady-state pressure drop across injector, AP’ 0. 2 

Nondimensional axial velocity difference, Av* 0.01 

Specific-heat ratio, y 1. 2 


The influence of these parameters on the calculated results can be found in reference 2. 


Chemical Reaction Model 

Stability limits for the chemical reaction-rate model are shown in figure 3 for a 
limited supply. With an unlimited supply all disturbances that were tested excited in- 
stability under all operating conditions values) because of the extreme sensitivity of 
the exponential function of the chemical reaction equation to a disturbance. However, 
when a realistic value for the concentration of unburned propellant is used (0. 1 or 0. 01 
mass fraction), the stability limits occur at very low burning rate parameters = 0. 1 
or less). As discussed in reference 2, these low burning rate parameters occur only 
for specialized small-scale research combustors. 


Vaporization Model 

Stability boundaries for the vaporization model are shown in figure 4. The boundary 
for infinite supply (M /u). t = a) is the same as that shown in reference 2 for a 
value of Av' of 0. 01. As the quantity of unvaporized material decreases, the stability 
limits are increased. (A larger amplitude of disturbance Pp is required to excite in- 
stability. ) Since 0 /u^ Q t is equivalent to the ratio of the time to burn an average 
particle of propellant to the wave time, it is seen that, when the burning time is much 
less than the wave time, for example, M y 0 A>j Q t = l/2rr, an extremely large distur- 
bance is required to excite the instability. 

The explanation for the increased stability with lower values of 0 Ax Q t w can 
be obtained from figure 5, which shows the time histories for the local burning rate, the 
unvaporized mass, and the pressure at two values of M v q /uk Q t w . With large values 
of 0 /aj. Q t the unvaporized mass remains constant during the disturbance, while 
for low values of Q /c 0 j Q t w the unvaporized mass decreases as the pressure and 
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burning rate increase. Since the burning rate as given by equation (12) is proportional 
to the square root of the unvaporized mass, the burning rate is reduced by the 
unvaporized-mass term during the first 1/4 cycle of the disturbance. This reduction 
results in less energy being added to the wave; hence, the disturbance decays. 


Atomization Model 

Stability boundaries for the atomization model are shown in figure 6. Stability 
lines are again shown for various values of M /o>. t . When the stability boundaries 
for the vaporization model (fig. 4) are compared with those for the atomization model 
(fig. 6), it is seen that in some regions the vaporization process is more sensitive to a 
disturbance and that in other regions the atomization process has the greater sensitivity. 
Therefore, both processes must be considered important in controlling combustion in- 
stability. 

Time histories of unatomized liquid, burning rate, and pressure are shown for the 
atomization model in figure 7. The histories are similar to those in figure 6 for the 
vaporization model; that is, for large values of M„ Vw, t the system is unstable, 
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Figure 5. - Time history of disturbance for vaporization model. Burning-rate 
parameter, 1; nondimensional peak-to-peak pressure disturbance, 0.01. 

and the unatomized liquid remains constant. For small values of M /ix>- t the 

a, o' l, o w 

system is stable, and the unatomized liquid decreases as the burning rate and pressure 
are increased. This decrease in unatomized liquid reduces the burning rate, which re- 
sults in insufficient energy to drive the disturbance. 


Vortex Flow Model 

Recent studies (ref. 3) have shown that with viscous forces a transverse pressure 
oscillation can generate a vortex within the combustor. Therefore to determine the 
significance of such a flow pattern, additional calculations were performed to determine 
the influence of a vortex on stability. These calculations were performed by using the 
vaporization model with a constant vortex velocity v yr added to the calculated values of 
gas velocity in the tangential direction, v^. Results of these calculations are illus- 
trated in figure 8. The stability limits are dramatically changed by a small vortex vel- 
ocity. For a given vortex velocity the system has the same stability at low burning rate 
parameters that it had for the vaporization model. At some critical burning rate param- 
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Figure 6. - Instability boundaries for atomization model. 


eters, dependent on the vortex velocity, the vortex model shows that almost any distur- 
bance will grow. Larger critical burning rate parameters are required for smaller vor- 
tex velocities. The effectiveness of baffles can be partially explained by the vortex 
model results. The vortex velocity can be eliminated by using baffles, and the stability 
of the vaporization model boundaries can thus be increased. 

The vortex model can also explain the results reported for the investigation of refer- 
ence 4, in which a tangential gas flow was used to excite instability. The quantity of 
tangential gas flow required to induce instability was dependent on the operating condi- 
tions of the combustor. 


EXPERIMENTAL RESULTS 

To test the validity of the theoretical results, an experimental program was also 
conducted.' Since the calculations were only performed for an annular segment of a com- 
bustor, an annular or torus combustor was used in the experimental program (fig. 9, 
p. 15). The torus had an 8-foot diameter and 10 injectors with nozzles equally spaced 
around the chamber. Four different types of injectors were used to change the operating 
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Figure 7. - Time history for atomization model. Burning-rate parameter, 1; non- 
dimensional peak-to-peak pressure disturbance, 0.02. 

conditions. The injectors (shown in fig. 9) were (1) an unlike-triplet injector with two 
hydrogen streams impinging on one liquid-oxygen jet, (2) an impinging-jet injector with 
two liquid-oxygen jets impinging on each other and four parallel jets of hydrogen, (3) a 
coarse parallel-jet injector with one jet of liquid oxygen and two jets of hydrogen, (4) a 
fine parallel-jet injector with four jets of liquid oxygen and four of hydrogen. Two-inch- 
diameter windows were placed in the combustor so that the combustion process could be 
observed and the steady combustion characteristics established. 

The torus combustor was operated over the following range of conditions: 


Steady-state chamber pressure, P £ , lb/sq in. abs 40 to 120 

Oxidizer-flow rate, lb/sec ’ 0. 95 to 3. 5 

Fuel-flow rate, lb/sec 0. 15 to 0. 45 
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Total flow rate, lb/sec . . 
Oxidizer -fuel weight ratio 
Oxidizer injection velocity, 
Hydrogen injection velocity, 


v i,ox’ £t / sec 
v H,,i’ ft/sec 


Oxidizer pressure drop, AP^ x 
Hydrogen temperature, °F . . 
Liquid-oxygen temperature, °F 


1. 2 to 4. 0 

2. 0 to 8. 0 
15 to 55 

300 to 600 

0. 2 to 0. 8 
32 to 80 
. . -320 


To excite instability a tangential blast tube 10 feet long with a 2-inch inside diameter 
was added to the combustor. A burst disk that ruptured at 1000 pounds per square inch 
was placed at the combustor end of the blast tube. The tube was loaded with 75 grams of 
ammonium perchlorate containing 25 percent by weight polybutadiene acrylic acid, which 
was ignited during the run and which provided sufficient gas to burst the diaphragm and 
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:er tube 
D. fuel tube 


Four 0.089-in. -I. D., 
equally spaced fuel tubes- 



25 

32 





Figure 10. - Experimental and theoretical stability boundaries for torus combustor. 

Burning-rate parameter, 0.015. 

produce a peak-to-peak pressure oscillation of 20 percent of the chamber pressure in the 
torus combustor. 

The experimental results are shown in figure 10, where the disturbance amplitude is 
plotted against the mass accumulation factor 27rM Q /w i Q t . For the unlike -triplet and 
impinging- jet injectors the 20-percent disturbance always damped. The instability of all 
runs with the two parallel-jet injectors indicated that the combustion noise produced dis- 
turbances large enough to excite instability. The noise level was assumed to be the same 
as that observed for the triplet and impinging-jet injectors during steady-state combus- 
tion (4 percent of chamber pressure), although larger disturbances were present during 
ignition and while flows were adjusting to their average rate. The instability was always 
initiated sometime after ignition but before average flow rates were achieved. 

The curves in figure 10 are cross plots of data presented in figures 4 and 6 for a 
burning rate parameter of 0. 015. This value of *£ is for the torus with a radius 
of 48 inches, an average burning rate per total length of combustor m of 0. 131, and a 
combustor contraction ratio of 409 (obtained by dividing the volume of the combustor 
by the product of the nozzle area and the distance from the injector to the nozzle). The 
theoretical curves show that the vaporization model has a slightly lower instability limit 
than the atomization model for the operating conditions of the torus combustor. The ex- 
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perimental value of the mass accumulation factor was obtained by calculating the amount 
of vaporized propellant as a function of length (ref. 5). These calculations agreed with 
the observation at the window that all of the propellant was vaporized for the impinging- 
jet and triplet injectors, while only a small fraction was vaporized with the parallel-jet 
injectors. Thus the parallel-jet-type injectors provided sufficient liquid in the combus- 
tor to significantly reduce the calculated disturbance amplitude required to excite in- 
stability. 

The experimental and calculated results shown in figure 10 indicate that good agree- 
ment exists between theory and experiment. Theory indicates that a disturbance greater 
than 40 percent is required to excite instability with the impinging-jet and triplet injec- 
tors, and experiments showed that a 20-percent disturbance would decay. Similarly, the 
parallel-jet-type injectors, which were always unstable, should require only a 2- to 
6-percent disturbance to excite instability. Operating conditions of the torus combustor 
did not completely agree with the conditions used in the calculations. However, these 
would not significantly change the calculated curves because the calculations are quite 
insensitive to Jf, v* , AP’, and y. The minimum disturbance is very sensitive to Av'; 
however, the value of 0. 01 for Av* is very reasonable in that it represents the turbu- 
lence level available as a driving source when gas and liquid velocities are very low. 

High-speed motion pictures taken of the sprays during instability indicated that the 
sprays were undergoing breakup due to crossflow at the high-pressure portion of the 
oscillations. It was impossible to determine the rate of breakup or vaporization, how- 
ever, because the liquid-oxygen sprays are not clearly outlined in the pictures; there- 
fore, the volume of liquid could not be measured accurately enough to obtain slopes and 
hence atomization or vaporization rates. The amount of liquid in the combustor varied 
during the cycle as predicted by the theory. 


SUMMARY OF RESULTS 

Theoretical calculations have shown that the physical processes of vaporization and 
atomization probably are the controlling parts of the combustion process in determining 
combustion instability. Chemical reaction is not likely to be important for most rocket 
combustors because there is insufficient mass or energy available in gaseous form to 
sustain an oscillation even though the reaction rate is very sensitive to pressure, veloc- 
ity, or temperature variations. 

Theory also indicates that the stability of a combustor is very dependent on the mass 
of liquid in the combustor and the existence of any vortex velocity, as well as the veloc- 
ity difference and the burning rate parameter (ref. 2). 



Limited experimental studies with a combustor similar to the one used in the theory 
have confirmed the stability regions. Experiments also confirm qualitatively the vari- 
ations that occur during an oscillation. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, February 9, 1965. 
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APPENDIX A 
SYMBOLS 


A c 

cross-sectional area of com- 
bustor, sq in. 

A d 

average surface area of drops, 
sq in. 

A i 

injection hole area, sq in. 

A t 

nozzle-throat area, sq in. 

yl 

combustor contraction ratio, 
A c /A t , dimensionless 

a 

average speed of sound in 
gases, in. /sec 

C D 

coefficient of discharge of in- 
jector, dimensionless 


average steady- state mass 
fraction of unburned gases, 
dimensionless 

c* 

characteristic exhaust veloc- 
ity, ft/sec 

D 

molecular diffusion coeffi- 
cient, sq in. /sec 

E 

activation energy, in. lb/Lb 
mole 

f(y) 

function of gamma, 

g 

acceleration due to gravity, 
386. 09 in. /sec^ 

fl 

viscous -dissipation parameter, 
'V*/ Rp c,oe 


k 

preexponential constant in 
Arrhenius equation, 
cu in. /(lb mole) (sec) 

& 

burning-rate parameter, 
Rm/jl 

M 

mass concentration, lb/cu in. 

M a 

mass concentration of unatom- 
ized propellant, lb/cu in. 

M 

r 

mass concentration of unre- 
acted propellant, lb/cu in. 

M 

V 

mass concentration of un vapor - 
ized propellant, lb/cu in. 

m 

average burning rate of propel- 
lant per length of combustor, 
fraction/in. 


molecular weight of propellant, 
lb mass/lb mole 

N d 

number of liquid drops per unit 
volume in combustor, 
l/cu in. 

P c 

chamber pressure, lb/sq in. 

P c, o 

steady-state chamber pressure 

c, max 

maximum chamber pressure 

P 

c, min 

minimum chamber pressure 

p p 

nondimensional peak-to-peak 
pressure disturbance, 

P „ - P 

c, max c, mm 
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P,p tank pressure, lb/sq in. 

P v vapor pressure of liquid, 

lb/sq in. 

AP* nondimensional steady-state 

pressure drop across injector, 



R radius of combustor, in. 

universal gas constant, 18 510 
(in. )(lb)/(°R)(lb mole) 

Re Reynolds number, dimension- 

less 

r^ average drop radius, in. 

r. average jet radius, in. 

Sc Schmidt number, dimensionless 

T temperature, °R 

t time, sec 

t„ breakup time of liquid, sec 

-D 

t w average wave time, 27rR/a, sec 

t* nondimensional time, t/t 

w 

U velocity difference between liquid 

and gas, in. /sec 

v gas velocity, in. /sec 

Av steady-state average axial ve- 

locity difference between 
liquid and gas, in. /sec 

v' nondimensional gas velocity, v/a 

Av* nondimensional axial velocity 

difference, Av/a 

v vr steady-state average vortex ve- 

locity, in. /sec 

v' nondimensional vortex velocity, 

v /a 
V vr / 


v axial gas velocity, in. /sec 

z 

Vg tangential gas velocity, in. /sec 

We Weber number, dimensionless 

z axial coordinate 

z’ nondimensional axial coordinate, 

z/R 

a correction factor for mass trans- 

fer, dimensionless 

y specific -heat ratio, dimensionless 

9 angular coordinate 

p gas viscosity, lb/(in. )(sec) 

p gas density, lb/cu in. 

a> burning rate, lb/(sec)(cu in. ) 

u> atomization rate, lb/(sec)(cu in. ) 

a 

ok injection rate, lb/(sec)(cu in. ) 

co chemical reaction rate, 

r 7 

lb/(sec)(cu in. ) 

u> v vaporization rate, lb/(sec)(cu in. ) 

Subscripts: 
f fuel 

l liquid 

o steady-state conditions 

ox oxidizer 

Superscript: 

' nondimensionalized parameter, 

equal to parameter divided by 
steady-state parameter, i. e. , 

T* = T/T q (except for velocities, 
which are divided by speed of 
sound, and time, which is 
divided by wave time) 
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APPENDIX B 


EQUATIONS FOR RATES AND MASSES 
Injection Rate 


The local instantaneous injection rate is given by 


w i “ C D A i V^T P c^ 2p i g 
The steady-state injection rate is 


"i, o = ‘lA ♦ / < P T- p c,o> 2 ^ 

«i- — -1 f- 

w. x p 
1, o f 


' P T - P c 


P™ - P 
T c, o 


r ap 


Atomization Rate 


The time to atomize a jet or drop of liquid in a cross current for values of 
We/Re l/ 2 much larger than 1 is (ref. 6) 

/o \ 2 / 2 / r 3/4 

t„ = 6.43^ U\ !i-V /4 


B 


w p, 


U (2We) 1,/2 


where Re = rj(Up/p.) and We = rj(pU /a). If the instantaneous atomization is the 
as the average for complete breakup, the atomization rate is 




u>a _ 


and 


(Bl) 


(B2) 


(B3) 


(B4) 

same 
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CO 


M 

a, o 


a,° t 


B, o 


then 


co„ 


“a = 


co 


= M> 
a 


a, o 


( p -) 5 / 12 ^- j 1 - 
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= M^(p’) 5/12 


v 2 + (Av ) 2 


(Av)‘ 


0. 625 


= M ’ (P*) 5/12 


.2 


n0. 625 


1 + 


vAv' 


(B5) 


Vaporization Rate 

The rates of vaporization for large Reynolds numbers as given in reference 1 are as 
follows: 


CO = 
V 


N d Dw f A d aP v 


^T2r d 


0.6ScV3[ 2 r d (v-v lo )]V2g) 


1/2 


(B6) 


oo. 


_ N d Dw £ A d,o aP v 


v ’° <^>T2r 


d, o 


0. 6 Sc 1/3 j2r dj 0 (v - v { , Q )J 


1/2 


(B7) 


oo 


00‘. = 


A /r X 1 / 2 / l/ 2 /v-v N 1 / 2 

A d / r d,o\ /p\ ( V 


- 


w v, o A d, o \ r d / v c 


lV o ' v £, o/ 


(B8) 


r, X 3 /2 / .1/2 / 2 + Av ^ 


1/4 


co* = 


d, Oy 




Av 
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= (M’) l/ 2 (p ') l/2 


1 + 


\Av* 


1/4 


(BIO) 


Chemical Reaction Rate 

The chemical reaction rate is given by a second-order rate expression in which an 
Arrhenius rate constant is assumed: 


CO 


M 

f xxA nir 

OX I 


— — = k eX p (dL ) 

m * m — m { 


(Bll) 


-E 


co M M t 

r,o = r, ox, o _r L f i o k exp 

m t m m f l^T 

ox, f ox f \ o/ 


(B12) 
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Mass of Unatomized Propellant 


The mass of unatomized propellant is determined by a mass balance: 


M = M + f (u>. - u) )dt 
a a, o Jq k i a' 


(B14) 




M' = 


M„ 


a, o 


CO. f ( CO. CO } 

= l + ^° / (- 1 

M a,o X \ W i,o w i,oy 


t dt’ 
w 


(B15) 


Assuming that all the propellant is consumed in the combustor gives 


CO = CO. = CO 

v, o 1,0 a, o 
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(B16) 


M £L 


= i + r (w? - co* )dt* 

M_ _ W v) a 


a, o 


Mass of Unvaporized Propellant 


The mass of unvaporized propellant is also determined by a mass balance: 


since 



CO ss CO. = Cl) 

v, o i, o a, o 



copdt’ 


Mass of Unreacted Propellant 


The amount of unreacted propellant is also determined by a mass balance: 


where 



- co r )dt 


M. 




p 

o*o 


(B17) 

(B18) 

(B19) 

(B20) 

(B21) 

(B22) 
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M* = l + -i-2 
r 


w. J- c o. 2irR 

° — = — = 2713? f(y) 


(See eqs. (A29) and (A31) of ref. 2.) 


‘r * 1 + f a V K - “r> dt ' 
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